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Rased on the complete Stokes wave model with the bias term at.d using a simple 
mapping approach and an iteration solution method, we established a formula for the 
joint probability density fuction of the surface slope-elevation of a nonlinear ran- 
dom wave field. The formula reguires three parameters to define the whole density 
function: the rms surface elevation and slope values and the significant slope. This 
model represents the dynamics of the wave in a more direct way than the Gram- Chari ier 
approximation. Based on this new statistical model and laboratory experiments, form- 
ula and numerical values of EM bias and dynamics bias are derived. The results inidi- 
cate that various biases should be considered seriously if accuracy of the altimeter 
measurement is required in centimeter range. 


1. INTRODUCTION 


Wind waves are always random. The randomness is the result of the generating 
forces as well as the consequence of the dynamic processes in wave evolution which 
induce different kinds of instabilities. Consequently, the only meaningful descrip- 
tions of the wind wave fields are the various statistical measures; among them the 
probability density functions of the surface elevation and slope are the most basic 
ones. With the recent development of active microwave remote sensing the techniques, 
the need for the statistical description of the ocean becomes more urgent, for the 
return signals of the radars are various convolution of radar signals and the ocean 
surface and these return signals are our only Information sources. Successful extrac- 
tion of the qeogphysical parameter, therefore, depends critically on our knowledge 
of the ocean surface statistical properties. For the remote sensing applications, 
the most important statistical measure is the joint, probability density function 
of the slope and elevation of the ocean surface. 


This specific probability density function Is necessary for evaluating the back- 
scattering equation for a near nadir-looking radar as stated by Rarrlck (1972), Brown 
(1978), Valenzuela (1978), and Jackson (1979). Unfortunately, due to the nonlinear 
nature of the ocean waves, a realistic functional form has been elusive. Faced with 
this difficulty, past Investigators have adopted a simple model that assumes all the 
statistical processes to be Gaussian and joint Gaussian. With this assumption, the 
statistical properties of the ocean surface become totally Independent of the control- 
ling dynamic processes. Consequently, the treatment of the ocean surface statistics 
Is no different from any other random Gaussian process such as the noise of electronic 
circuitry. Indeed, most of the p r esent statistical results used in describing the 
ocean surface can be traced to the classical set of papers by S. 0. Rice (1944, 1945 
and 1948), with certain extensions Into two dimensions by Longuet-HIgglns (1957, 1963, 
for example). Because the randomness of the ocean surface wave Is the consequence 
of the generating forces as well as the controlling nonlinear dynamic processes which 
Include various types of Instabilities (see, for example, Phillips, 1977 and Yuen 
and Lake, 1982), the resulting surface geometry can not be adequately described by a 
linear superposition of statistically Independent events. Thus the statistical pro- 
cess can not be Gaussian. Therefore, any results based on a Guasslan assumption. 
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though perhaps useful for some qualitative purposes, do not satisfy the need of our 
present and future engineering and remote sensing requirements. 

The nonlienar effect gravity waves on statistical properties of the wave field 
has been investigated first by Longuet-Higgins (1963), who treated the surface statis- 
tics as nearly Gaussian and represented it by a Gram- Chari ier series ’n which the 
skewness nf the surface elevation was determined through rigorous dynamic equations. 
Thi' analysis was later extended by Jackson (1979) to the joint probability density 
distribution of slope and elevation. These results represented a major breakthough 
in the non-Gaussian statistical description of the ocean surface. Successful as they 
were, there are certain shortcomings: the approximation gives negative density values 
and it requires higher moments to implement. These shortcomings put a limitation the 
applications of the results. 

Recently we established a mapping technique to model the non-Gaussian process of 
the ocean surface. This technique was first used by Tayfun (1980), and was later suc- 
cessfully modified and applied to the non-Gaussian processes described by the nonlin- 
ear surface wave elevation density function (Huanq et aU 1933) and joint slope and 
elevation density function (Huang et aj_. 1984). In this paper we will summarize the 
technique used in constructing the non-Gaussian density functions and discuss the 
specific applications to the remote sensing techniques. 


?.. THE NEW STATISTICAL MOO EL 

The new statistical model for a nonlinear deep water wave field is based on the 
Stokes expansion, i.e. 

( = ak + a cos % -*• &>sz% + - cx>s30( + — . (2.1) 

Let us start by considering the linear Gaussin case for which the sea surface can be 
represented by 

C 1 = a cos x , (2.2) 
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with x as the phase function 

x = kx - nt | 

where k is the wave number, n is the frequency and 0 is an arbirary phase shift uni- ! 

formly distributed. It is well known that if a is Rayleight, the cl is zero-mean j 

Gaussian (see, for example, Papoulis, 1965). For the same amplitude and phase func- ; 

tion, | 

c2 = a sin x , (2.3) 

Is also Gaussian. Furthermore, r 1 ad c 2 are statistically orthogonal to each other. 

Basrfd on c 1 and r, ?., we can define the normalized random variables 

Z, - fz . Zz _ a si (2.4) j 

with the overbar Indicating mean quantity, and the joint probability of l\ and Z£ Is 

simply ! 
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(2.5) 

If the surface wave profile is represented by the Stokes wave, a normalized zero- 
mean random variable can be defined as 

„ . -4z I = ,3-? 

( ( £2 _ (T ’ (2.6) 

where a is the standard deviation of the surface elevation. From (2.1) we can also 
get the surface slope 

q 

— - ak S»VlO( - <k z k z smZX ~ y * 3 J< 3 Sm 3% - ... . (2.7) 

Aqain, we can define another normalized zero-mean random variable from (2.7) as 

\ ~ $*/ ( Sjf V /z (2.8) 

In terms of Zi and Z£, n and r become 

s i - - 3 - ^ - i 

in which 

N = ( I-+ «Tt 2 § 2 )Yz , 

M - ( l-Hfe-Tt 4 f 

and § is the significant slope of the wave field defined as 

§ - l!‘)V» , 

wl th X =2n /k . 

Since we know the joint density function of Zi and Z 2 , it is easy to write the 

kSlwJ -! n ?i ty / U 2 Ct1 °I ? f ^ ny tW ° I unCti ' * of Z i* and Z 2 through a mapping technique 
known as the fundamental theorem of probability fPapoulis, 1965), i.e. H 

f(z,,z 1 )d2,j2j=- p {z,( 1 . 5 ),z i < rt )j 3 (i^ )dl<js ( (2.10) 

in which J() Is the Jacobian of the transformation. 

densItTfunctlon^s:^^’ — — (1984) obta1ned the slope-elevation 

“ ~{l-4 n <rK4 [(«qVz <JM z $VjN*]<r 4 k 2 } 

r N rj 4- M 2, + 4 r| m 2 3 o-lc “l 
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dl strl but lo^o / 16 j0lnt specu1ar po1nt and elevat * on distribution Is the conditional 


For a long crested wave field or a narrow band of energy spread case. 


( 2 . 12 ) 
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Thus a good approximation of the specular point distribution can be obtained by simply 
setting r =0 in (2.11). Hence, in terms of the normalized variable. 


V I S*° ) = [l'4^<rk + i~ *] 2 - 4 ) 3-^k 2 j 


• exp 



N 2 »| z N z <rk 'l ) 


( 2 . 14 ) 


where a is the normalized factor to guarantee 

CD 

! ?(<\ i % *o ) ^ = 1 • 
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From equation (2.11) and 2.14) we can see that the density functions are obvious- 
ly skewed with respect to n. This skewness is caused by the harmonic distortion of 
the wave profile described by Stokes expansion. Comparison between the throretlcal 
and experimental results have been presented by Huang et£l_. (1983, 1984). The agree- 
ments are all good. 


3. REMOTE SENSING APPLICATIONS 


A particularly Interesting aspect of the specular point density functions is the 
location of the mean. More specifically, the difference between the mean of the spec- 
ular point density and the mean of the surface elevation density functions. This 
difference Is the critical electromagnetic bias (Jackson, 1979; Walsh et al., 1983), 
for the mean rea surface measured by any nadir looking electromagnetlc^eTTce Is not 
the true mean out the mean of the reflecting facets or the specular points. Based 
on our observational data, the mean Is determined and plotted as a function of the 
significant slope In Figure 3.1. The data show an almost linear Increasing of the 
bias. The mean value of the theoretical model can be calculated by the ratio of the 
first moment to the zeroth moment of the density function. The value of the bias Is 
the difference between the corresponding values of the mean In specular point density 
and elevation density functions. The numerical value calculated based on (2.14) Is 
also given In figure 3.1. Although they do not exactly conlnclde, the trend Is 
unmlstakeable. Better model with more realistic wave profile In two-dimension may 
be needed. 


Another aspect of the nonlinear effects In the wave profile Is the non -zero 
mean - the deviation of the local mean surface under waves from the undisturbed calm 
surface. The existence of this non-zero mean In a nonlinear wave train h * been shown 
by Stokes (1880), Rayleigh (1917), De (1955) and, more recently, Schwartz (1974) and 
Longuet- Higgins (1975). In the past, the terr a"k/2 In equation (2.1) has been treated 
as a negligible constant In all wave studies. It Is discarded by simply shlftlnq 
the coordinate system up by the exact amount of a’ k/2. This Is permissible as long 
as our Interest In the oscillatory part of the motion. But the coordinate shifting 
Is no longer permissible for our present interest because the amount shifted Is 
precisely the dynamic sea state bias which could have far reaching consequence In 
the accurate determination of the mean sea surface using an altimeter. 

The existence of this non-zero mean in a nonlinear wave train can be explained 
by mass balance in the wave motion. It has been shown by Stokes that a non-zero 
mean mass transport always accompanies the wave motion. For deep water waves, this 
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FIGURE 3.1 EH Bias as function of significant slope. 

— numerical value according to the new statistical model 
0 Laboratory simulated data. 




velocity is 


UjL = ofo ^ 


(3.1) 


where Z is the vertical position of a particle. The rate of total mass transport 
through any vertical section is given by 


Q = { an k e zkz dz = an /z 


( 3 . 2 ) 


The surface elevation change caused by this mass transport can be calculated as 


&%. v , 


Q T/ A 


(3.:) 


in which T and X are the wave period and wave length respectively. Using the defini- 
tion of wave number and frequency, one can easily show the equivalence between (3.3) 
and t 'e constant term in (2.1). 


For a random wave field, the Stokes drift Is given by Huang (1971) as 

a, * S f 2nk 

fern ~ 

Therefore, the rate of the total mass transport is 
Q * \ ° u x dz 
and the bias is '*** 

n _ \° [ f !~! 2mk tftV.o) dz . 

= j f k X(K.m) dk dn , 

k a 


(3.4) 


(3.S) 


(3.6) 


in which X(k,n) Is • * <:t1onal wave number frequency spectrum. In order to evalu- 
ate (3.6), a sped >< s r -urn form has to be adopted. Since the quantity sought here 
In a scaler, Ah, we can reduce (3.6) by using the dispersion relationship to eliminate 
k and reduce X(k,n) to a simple frequency spectrum. Then adopting the simplified 
Wallops spectrum model (Huang et al., 1981), we obtain 


* mi § (m-.',/( W . }) , (3.7) 

where m is a cimple function of the significant slope, §, defined by 

m - | I , 

and 

with x n as the wave length of the wave having a frequency at the peak of the spectrum, 
n 0 . 

If tne significant wave helqht, Hj/ 3 , is used as In most practical applications, (3.7) 
can be written as 




(3.8) 




The typical value of § in the open ocean Is around 0.01, which qives the value of m 
around 10; the dynamic sea state bias, Ah, would be ground 1.5% of the significant 
wave height. During a storm, the § value could reach 0.02 as an upper bound, and then 
Ah would be as high as 3% of H1/3. For the sea state dominated by swell, however, 
the dynamic bias can almost be neglected because § would be of the order of 0.001 or 
less. 


As an estimate of the magnitude of this bias, we use the empirical formula qiven 
by Neumann and Pierson (1965), then 

H 1/3 = 2. 12xl0- 2 W 2 (3.9) 

for a fully developed sea. Under such conditior.s the phase velocity of the energy 
containing waves equals the wind speed. Thus the sea state bias can be written as 

= 2. IS *IO~ A W 2 ( m-0/(rv)-n . (3.10) 

This value is plotted in Figure 3.2 For all wind speed, the high sea s^ate this 
bias is not negligible. 

4. CONCLUSION 

In this study, we summarized the non-Gaussian statistial models based on Stokes 
Wave expansion. The rew statistical model indicates the existence of skewness which 
is the cause 0* the EM bias in the remote sensing application. Furthermore, the exis- 
tence of the non-zero mean in Stokes expansion which contributed to the Skewness can 
also cause a dynamic bias. The magnitudes of both bias terms are estimated. For 

future altimeter applications, we should seriously consider the correction of these 
biases. 
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